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Abstract 

We propose a field theory model for dark energy and dark matter in interaction. Comparing 
the classical solutions of the field equations with the observations of the CMB shift parameter, 
BAO, lookback time and Gold supernovae sample, we observe a possible interaction between dark 
sectors with energy decay from dark energy into dark matter. The observed interaction provides 
an alleviation to the coincidence problem. 
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Recently there have been several papers dealing with interacting Dark Energy and Dark 
Matter [l, 2, iQ,^. It was argued that Dark Energy and Dark Matter interact via a small 
coupling, of the order of magnitude of the fine structure constant [6]. Employing some data 
sets from observational cosmology, including CMB shift parameter, BAO, age parameter 
and supernovae observations etc., it has been shown that the interacting model is a useful 
and robust model at the order of one standard deviation , while some observations 

are good enough at two a level, providing some confidence on the results. 

Most available discussions on the interaction between dark sectors are concentrated on the 
phenomenological investigations. It is of great interest to describe the interaction between 
dark energy and dark matter from a fundamental field theory point of view. Recently, some 



attempts have been proposed in 



13j. In order to follow this thread, we consider now an 



interacting field theory with two fields describing each of the dark components, a fermionic 
field for Dark Matter and a bosonic field for the Dark Energy, which here we adopt to be 
the tachyon field sl- 11|. We thus consider the Lagrangian 



C 



where a is a constant with dimension Me^~^, (3 a coupling between dark energy and dark 
matter fields, V{(p) the tachyonic potential and g the determinant of the metric. For a 
Friedmann-Robertson- Walker cosmology g^u = diag{l, —a{tY, —a(t)'^, —a{tY) one finds the 
equation of motion for the scalar field to be 
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-. We also have 

a 



dt 
di 



0, 
. 



(3) 
(4) 



From the latter, \I'\I' 



. We note that such a result follows from the homogeneity 



assumed for Dark Matter distribution. Thus, Dark Matter in our model just follows the 
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universe expansion, what is consistent with the cosmological principle. Moreover, 
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where we defined the effective mass M* = M — f3(p. Note that uJ^p 
Deriving [5] and [7] with respect to time and using [2] and HI we get 
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These equations are very similar to those usually used as a phenomenological model for the 
interaction between dark matter and dark energy 0, Q, 12]. The right hand side in the above 



equations does not contain the Hubble parameter H explicitly, but it does contain the time 
derivative of the scalar field, which should behave as the inverse of the cosmological time, 
replacing thus the Hubble parameter in the phenomenological models. 
The Friedmann equation for a fiat universe reads 



where M^^ 



Some analytic solutions in the pure bosonic case have been found in 
potential 
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and 
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for the 



(12) 



We choose, at this moment, n = 2, which leads to a power law expansion of the universe. 
However, we shall see that this choice is really not important and some properties depend 
little on the actual choice of n. This actually lowers the appeal of the present model. 

Let us now compare the interacting DE with the observational data. We will compare the 
interacting tachyonic model with the luminosity distance of the Gold supernova sample (182 
type la supernovae observations), the shift parameter of CMB radiation, the measurement of 
the Baryonic Acoustic Oscillations (BAO) and ages of galaxy clusters (see j3, 3, 14, 1^, 16|). 



It is convenient to rewrite ([2] 
with n = 2, ([2]) becomes 

^ijJ\Juj + 1 



in terms of two first order equations. Using V^(v?) as in (fT2 
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QHioiio + 11 
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Above, we defined = Equating the actual values of and with the observed 
values, = SMp^ifQ^l^Q and p^^^ = SMpiH^^l — Q^^), we can replace 'i'o^o and m by ob- 
servable quantities: M^^o^o = SM^^ and am^ = 3M^iH^{l-n^^)(pl^/\^\. 
Therefore, 
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where Hq 
and 



2.133/1 X 10 "^^MeV is the value of the Hubble parameter today, fi^^ 
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The parameter a is fixed as -v/a = 1.607 x 10^^Me\^ ^ ~ [Hn x MeF) ^. Since Sf 

n n 

|9| [lOl, such a was chosen such that the last term in (fT5|) is of the same order of magnitude 
as the other terms. In fact, we shall see - and it is easy to infer from the above equations 
— that only can be obtained, that is, a and M can be absorbed in the redefinition of 
j3. Thus we have, as parameters of the model, ( jv^-^ , 0O; h,VL^,^^,ujQ). 



In 



171], the lookback time method was discussed, allowing to use the cluster age to fix 



the parameters. Given an object i at redshift 2,, its age t{z.i) is defined as the difference 
between the age of the universe at Zi and the age of the universe at the formation redshift 
of the object, zp, that is. 



dz' 



dz' 



where ti is the lookback time given by 

Using ( fT9l) . the observational lookback time t°^'^{zi) is 

= tL{zF) - tiz,) = [C - tiz,)] - iC - tLizp)] 

= C - t{z^) - df , (21) 
where t'^''^ is the estimated age of the universe today and df is the delay factor, 

df = C-tL{zF) . (22) 

We now minimize Xaty 

i=l + 

where tL{zi,p) is the theoretical value of the lookback time in Zi, p denotes the theoretical 
parameters, t°^'^{zi) is the corresponding observational value given by (!2T!) . ai is the uncer- 
tainty in the estimated age t{zi) of the object at Zj, which appears in fl2T]) and a^oha is the 
uncertainty in getting t'^'^ . The delay factor df appears because of our ignorance about 
the redshift formation zp of the object and has to be adjusted. Note, however, that the 
theoretical lookback time does not depend on this parameter, and we can marginalize over 
it. 

In [l^ and the ages of 35 and 32 red galaxies are respectively given. For the age 



of the universe one can adopt t°^' = 13.73 ± O.UGyr [20|. Although this estimate for t°'" 
has been obtained assuming a ACDM universe, it does not introduce systematical errors in 
our calculation: any systematical error eventually introduced here would be compensated 
by the adjust of df, in (l2Ti) . On the other hand, this estimate is in perfect agreement 
with other estimates, which are independent of the cosmological model, as for example 



= 12.6j;2;4Gyr, obtained from globular cluster ages [2l| andtg = 12.5±3.0Gyr, obtained 



from radioisotopes studies 



'221. 



For the cosmic radiation shift parameter in the fiat universe we have 



where zis = 1089 is the last scattering surface redshift parameter. The value R has been 



estimated in 



23| from the 3-years WMAP [2^ results as Robs = 1-70 ± 0.03, for the flat 



universe, and is very weakly model dependent. Thus 

2 [R — Robs]'^ 
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Baryonic Acoustic Oscilations (BAO) [25| is described in terms of the parameter 

1 rZBAO rl'y' ^ /o 

zbao Jo E{z')' 

where zbao = 0.35. It has been estimated that Aobs = 0.469(7^)-°-35±0.017, with = 0.95 



24| being the scalar spectral index. We thus add to the term 



XU = ^^^^^ . (27) 



Finally, we add the 182 supernovae data from SNLS [26|, recent supernovae from 
HST/GOODS and further old data, as compiled by Riess, et. al. 0]. Defining the distance 
modulus 

fi{z) = hlogio [c(l + z) / — -] + 25 - Uog^oH, , (28) 



we have the contribution 



182 



We use que expression 
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.=1 



= xIn + xIao + xImb + xL + + , (30) 



where the last two terms correspond to gaussian priors for h [28j] and Qqi^ 27|], respectively: 
hobs = 0.72 ± 0.08 and Qnobs = 0.28 ± 0.04. 
The likelihood function is given by 

r( P ^ .o ^ \ xHjJ^,^o,h,n^,,ujo) 

^(tT^' '^0' ^*o) ^o) oc exp[ . (31) 

We present in the table the individual (marginalized) best fit for each parameter, with 
respective deviations. Figure (la) shows the curve fi{z), corresponding to the global best 
fit. Figure (lb) shows the fit of the lookback time tL{z). 



Table 1: Values of the model parameters from lookback time, BAO, CMB and SNe la. 



h 



0.631 ±007 ±0.015 ±0.022 



*0 



0.324 ±0.017 ±0.033 ±0.050 



UJQ 



_n Q7Q+0-106+0.227+0.283 
U. J ( -"-0.003-0.013-0.018 



boHo > 1.3 (2a) 







ji^km, < 0-071 (2a) 
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FIG. 1: Results corresponding to the global best fit (lookback time + CMB + BAO + SNe la). a) 
Theoretical distance modulus compared to 182 SNe la data, b) Theoretical tiiz) versus 67 galaxy 
clusters data. 



The parameters 0o and are strongly degenerated. Indeed, in figure 2 we see that the 
effect of both parameters are the same on the densities of Dark Energy and Dark Matter: 
increasing 0o is equivalent to decreasing j^j^- Here, it would be convenient to observe that, 
even in the non interacting case (when = 0), as we can see in figure 2b, the fi^, never 
goes to one. Rather, the ratio of the Dark Matter to Dark Energy densities remain constant 
in the Dark Matter domination era, because, in this era, the equation of state parameter 
of the Dark Energy approaches zero, and thus Dark Energy behaves as Dark Matter in 
this period - this feature had already been underlined in [lO|. Therefore, in the Tachyonic 
Dark Energy model the coincidence problem is less serious. However, the ratio of the Dark 
Matter to Dark Energy densities depends on the parameters of the model. In particular, 
lower values of 0o turns this ratio higher. The introduction of the coupling furnishes an 
additional improvement in the coincidence problem, diminishing the Dark Matter to Dark 



7 



Energy ratio, if the coupling constant is negative, as we can see in figure 2a. Positive 
values of jf^ -, on the other hand, aggravates such a problem. 





In(a) 
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FIG. 2: Relative densities of Dark Energy and Dark Matter, Jl^ and as functions of the scale 
factor a. a) for constant {(PqHq = 2.5). The dot-dot-dashed, dot-dashed, solid, dashed and 



dotted lines are for 

/3 



+0.125, 



-F0.0625, 

/9 



/3 



0, 



-0.125 and 



My^aHo ~ ~0.25, respectively, b) for constant {j^^ = 0). The dot-dot-dashed, dot-dashed, 

solid, dashed and dotted lines are for (poHo = 1-8, (PqHq = 2.0, (poHo = 2.5, (poHo = 5.0 and 
= 7.5, respectively. 

In order to compare the model with our previous predictions we compute 

the likelihood functions concerning the various parameters of our model. Our main previous 
prediction concerns the behaviour of the interaction, especially its sign. As it turns out, the 
model is very degenerated, but most of the allowed values of j3 are consistent with a negative 
coupling. 

In figure 3 we plot the behaviour of the f3 versus 0o contour for la and 2a. Since the 
diagram is unbounded for negative /3, its allowed values are, generally speaking, negative, 
although we cannot rule out a small positive coupling. The likelihood of /3 marginalizing 
all the other parameters is not normalizable, being consistent with all values of j^j^ below 
a small positive value. Most of the allowed values are negative, see figure 4. These results 
indicate that if there is a coupling connecting the dark sectors, it is more probable for the 
dark energy to decay into dark matter, which is consistent with the fact obtained in the 
study of thermodynamics 29|. 

In figure 5, the first diagram concerns the Dark Matter fraction versus j^j^ contours. 
We see that the observed value holds almost independent of j^y^ if this latter is not large 
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positive. Similar conclusions can be drawn for the diagram of the Hubble constant compared 
to 00 as well as Dark Matter versus (pQ. The likelihood of (pQ is shown in figure 6. 




-1,0 -0,5 

p/(Moc°=H„) 



FIG. 3: Two dimensional distribution of f3 and <j)Q {la and 2a contours). Notice that there 
is a strong degeneracy. P can go to arbitrarily large negative values. For positive values the 
function decays quickly to zero. The expectation of /? cannot be computed due to the fact that 
the distribution does not approach zero. Thus, /? should most probably be negative. 
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FIG. 4: The f3 likelihood function shows a behavior confirming the speculations arising in the 
previous figure. It is most probable that it is negative. 



In order to further understand this problem we consider the consequence of our formula- 
ion for an effective fluid interaction such as the one considered in several previous papers 
. In those works, the interaction term in the fluid conservation equations is of the 
form 6Hp, where S is the coupling constant. After some simple manipulations and using the 
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FIG. 5: Two dimensional curves displaying the probability distributions of (3 versus ri^Q, (po versus 
h and (po versus ^l^g, respectively. 
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FIG. 6: The (po likelihood function. 



definition of w, the equations iQ and ffTOl) can be put in the form 

/3 



p^ + 3Hp^{uj + l) = ^ VJTTpv^o(l + ^)' , (32) 



1 - 



/3 .,, 
/3 



p^ + 3Hp^ = ^ y^7TTp*o(l + , (33) 

where p^t^ = SMpiH^fl^^^. As we mentioned above, ^/a ~ H^^. Moreover, for 2; > 0, 
uj rapidly approaches very small values, so the interaction term in the r.h.s. of fl32|) and 

(!33|) will be of order of ( ]ifoP^o(l + z)^, very similar to the phenomenological 

interaction. Therefore, the phenomenological coupling constant 6 would be constituted of 

/3 



two theoretical parameters: and 



On the other hand, notice that, in the model considered here, the parameter 5 ^ g , 
is in fact an effective coupling constant. This appears in fl32l) and fl33|) and in the last term 
of (HM . The only other place where j^y^ appears is in (ITS]) , but in fact, H is not much 
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affected by tlie coupling, since decreases rapidly to a certain value (typically about O.80o), 
thereafter remaining constant (remember that we are integrating the equations "backward"). 



Therefore, 




FIG. 7: The Ukelihood function (IMD . 




6/H„ 

FIG. 8: Two dimensional curves displaying the probability distribution of 5 versus (j)^. 




I -| 1 1 1 1 r 1 T 1 1 1 1 

-1,00 -0,95 -0,90 -0-85 -0,80 -0,75 -0,70 



FIG. 9: Two dimensional curve showing the behaviour of (/>o and w. 
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FIG. 10: Dependence of the model on the value of n. 



We compute now the likelihood of the function 6 = — . The likelihood of S is 



determined from the likelihoods of j^y^ and (po accordingly to 

PM = -*)P(^,^„M(^M^„ , (34) 

^^^^^ ^(m7^'<^o) ~ P^^or{j^,(f)o)C{j^,(f)o) for some prior we choose for and 0o- 
We find the result shown in figure 7. It is convenient to introduce some remarks. 

As previously mentioned, we had rewritten ^'o^o in terms of observable quantities, 
\l/o\I/o = jvf (i''' ^ ) ' ^^^""^ ^0^0 is the number density of particles of Dark Matter, it 
must be positive. Thus we must have j£^4>o < 1- This constraint has to be included in 
the prior on 0o and jj^- Moreover, as a consequence, the effective coupling 6 has the same 
sign as f3 and more negative (positive) values of S are equivalent to more negative (positive) 
values of (3, turning the ratio of the Dark Matter to Dark Energy lower (higher) in the Dark 
Matter domination era. Thus, 6 has the same degeneracy with 0o as f3. However, S is much 
more restricted than j3, being contained in a small interval, as shown in figures 7 and 8. The 
constraint ]yf^ 4>o < 1 is equivalent to 0o < ~|) for 5 < 0, as it can be seen if we eliminate 
j^f^ in favor of 6 in the expression for 'i'o^E'o- Such a constraint can be clearly seen in figure 
8. The maximum of P{S), figure 7, depends on the prior on and (heavily) on the prior 
on 00, but the integrated probability for 5 < is very weakly prior dependent. We have for 
the probability for negative coupling P{6 < 0) = 90%. Thus our result is consistent with a 
small negative value of 6, which implies in DE decaying into DM, alleviating the coincidence 
problem. 

We also learn, in this connection, that in this model a more negative equation of state 
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for Dark Energy is connected with a larger value of 0o as shown in figure 9. However, such 
a result is rather model dependent. 

The first attempt of this study was to describe a few details about the Dark Energy and 
Dark Matter behaviours and thus we chose the model by specifying the index n to be 2. 
However, in figure 10 we see that the numerical solution shows nothing extraordinary for 
such a choice of n. Indeed, for n = 1, 2, 3 the results of the likelihood of the coupling (3 are 
surprisingly similar. 

All n tested in this work (until n = 10) were capable to reproduce the features of the 
observed universe (at least, of the background), and, in fact, the solutions of the equation of 
motion for all n tested had the same qualitative behavior: it reproduced the actual period of 
accelerated expansion, driven by the Dark Energy domination, and the Dark Energy equation 
of state parameter approaching zero in the Dark Matter domination era, forcing the ratio 



P^/P(f> to be a constant in this era. In [ll| two approximate solutions had been found. 



valid for a tachyon dominated universe: one which corresponds to ^ — 1, for < n < 2 
and another which corresponds to ~ 0, for n > 2. The numerical solutions for the exact 
equation of motion, encountered by us, in fact reproduce these predicted behaviour, but only 
asymptotically, in the far future (a(t) >> 1), when Vt^ ^ 1. These results are consistent 
with those found in ref jl3| , namely a dynamical attractor behaviour for general values of n. 

The model has further features that we consider being drawbacks, as e.g. the fact that 
w > —1 oi also the extreme non linearity of the action, rendering the calculation clumsy 
and the particle interpretation unclear. From a positive side, the problem can be opened up 
for more realistic models of Quantum Field Theory, as in 30 1. 

In spite of the simplicity of the model, the comparison of the model with the values of 
cosmological parameters leads to the conclusion that the interaction is consistent with the 
observations at least at one standard deviation, possibly at two, for negative coupling. This 
encourages us to look for more sophisticated theoretical field models for the explanation 
of the Dark Matter and Dark Energy behaviors as well as their origins. In particular, the 
present model does not account for equations of state with uo < —1. We could try to mimic 
this fact by taking more sophisticated potentials (V < 0) or a larger number of fields, but 
that would enlarge the number of parameters of the model. Since a transition redshift is 
not yet well established we prefer to stay with the more conservative case. 

A further point we have not dealt with is the comparison with the structure formation. 
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The procedure might imply further (and better) constraints for (3 and 0o- However, we did 
not find it worthwhile pursuing further this simple model and we left aside this possibility. 

As a conclusion we can state that it is reasonable to expect that DE and DM interact via 
a small but calculable and observable coupling, possibly giving an alternative to the usual 
cosmological constant explanation of Dark Energy. 
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